#include "Matrix.h"

using namespace std;

class ActiveRegion
{
  private:
//  number of grid points
    int grid_size;
//  holds the x (position) values of the grid points
    double * x_grid;
//  spacing between grid points
    double delta_x;
//  potential at each grid point
    double * V;
  public:
    ActiveRegion(char *);
    ~ActiveRegion();
	void MakeHamiltonian();
	Matrix Hamilt;
	double* Fork0;
}; 
